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Pion absorption and rescattering in the ANP model revisited 

I. Schienbein a and J.-Y. Yu b 

a Theoret. Physik II, Univ. Hamburg, Lumper Chaussee 149, 22761 Hamburg, Germany 
b Theoret. Physik III, Univ. Dortmund, 44221 Dortmund, Germany 

Single pion leptoproduction in the region of the (3, 3) resonance is currently of high interest for at least two 
reasons: (i) These reactions constitute an important part of the total cross section in low energy reactions and 
are utilized to detect neutrino oscillations in current and future long baseline experiments, (ii) Intranuclear 
rescattering of the pions in heavy nuclei results in interesting and sizable modifications of the free nucleon cross 
sections which are testable in electroproduction experiments. In this article we give a basic introduction to the 
pion multiple scattering model of Adler, Nussinov, and Paschos (ANP) with special emphasis on pion absorption. 
We also estimate the probability of multiple scattering. 



1. Introduction 

Neutrino nucleus scattering will be an impor- 
tant tool in future long baseline (LBL) experi- 
ments to precisely measure neutrino parameters 
like neutrino masses and mixing angles. These 
experiments will use neutrino beams with ener- 
gies of the order 0(1 GeV) where quasi-elastic 
reactions (QE), single pion resonance production 
(RES) and deep inelastic scattering (DIS) are all 
important. Moreover, due to the small neutrino 
cross sections heavy targets like oxygen, argon or 
iron have to be used. 

In the following we will deal with single pion 
resonance production, i.e., with the reactions 



T ->• Z + T' 



r±.° 



(1) 

Ar 40 , 26 Fe 56 ) 



where T is a nuclear target (sO 16 , 
and T 1 a final nuclear state. 

At this stage we make the basic assumption 
that the reaction Q can be described by two in- 
dependent steps pQ: 

1. Single pion production in vN scattering: 
In this first step a pion is produced by 
the scattering of the incoming neutrino off 
a nucleon in the target, such that nuclear 
corrections due to the Pauli principle and 
the Fermi motion of the nucleons should be 
taken into account. 



2. Multiple scattering of pions: 

Once the pion has been produced it will 
travel through the nucleus. During this 
journey the pion can have several reseat - 
terings in which it can change its charge or 
be absorbed. 
Step 2 can be described by a 3 x 3 charge exchange 
matrix M which depends only on the properties 
of the target [modeled by a charge density profile 
p(r)\ but which is independent of the identities of 
the leptons in step 1. It should be noted that the 
above assumption of two independent steps gen- 
erates predictive power, since the formalism can 
be applied to charged current (CC) and neutral 
current (NC) neutrino- and electroproduction of 
pions in the resonance region. 

Absorbing the Pauli suppression factor of step 
1 into the normalization of M, the measur- 
able final distributions (cross sections) of pions 



(7T+ 

butions (7r + 
pie form 

/ 



)f 



can be related to the initial distri- 
7r , ir~)i for a free target in the sim- 



da( ; 



T A :n+ 



) \ 



dQ 2 dI^ 
d<j( z T A ;ir°) 

&Q 2 dW 



\ AQ 2 dW ) 

V 

nuclear target 




(2) 



free nucleon 
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with 

da(JV T ;±0) = da(p; ±0) da(n;±0) 
dQ 2 dW dQ 2 dW { ' dQ 2 dW 

where the free nucleoli cross sections are aver- 
aged over the Fermi momentum of the nucle- 
ons. The matrix M depends on the target ma- 
terial and the final state kinematic variables, i.e. 
M = M[T;Q 2 ,W}. 

In the following we will restrict the discussion 
to isoscalar targets and refer the reader to ref. 
[2] for targets with a neutron excess. The charge 
exchange matrix M for isoscalar targets can be 
parameterized by three independent parameters 
Ap, d, and c: 

(1 — c — d d c \ 

d I -2d d . (3) 

c d 1 — c — d I 

So far there is no direct experimental informa- 
tion on these three parameters available. On the 
other hand a model for pion multiple scattering 
by Adler, Nussinov and Paschos (ANP model) 1 
allows to calculate these parameters and (aver- 
aged) charge exchange matrices for oxygen, argon 
and iron targets can be found in 

These matrices are useful to calculate CC and 
NC single pion production in the (3,3) resonance 
region induced by neutrinos j.'il lloKij and elec- 
trons and the rescattering effects are predicted 
to be quite large. For example, the pion energy 
spectra for w° 's produced on oxygen targets in NC 
neutrino production will be reduced by up to 40% 
whereas the analogous spectra for the charged pi- 
ons remain almost unchanged compared to the 
free nucleon cross sections 0. Similar results are 
predicted for the VF-distributions for 1-pion res- 
onance production in electron-oxygen scattering 
. For the 7r° the VF-distribution will be reduced 
by about 50% while the production of 7r + remains 
almost unchanged compared to the free case since 
the decrease due to absorption is widely compen- 
sated by an increase due to charge exchange. This 
marked signal of rescattering effects should be 
easily testable in low energy electron scattering 
experiments, for example at JLAB. 

The rest of this paper is organized as follows. 
In the next section we give an introduction to the 



ANP model. In Sec. El we discuss the normaliza- 
tion of pion absorption which we fix by experi- 
mental information. In Sec. 01 we discuss the dy- 
namics of the ANP model in more detail. Finally 
in Sec.|S]we summarize the main results. 

2. The ANP model 

In this section we introduce the multiple scat- 
tering model by Adler, Nussinov and Paschos P^. 
Due to space limitations the presentation will be 
short and we refer the reader to Ref. j: for a 
thorough discussion of the omitted details. 

The main ingredients of the model are as fol- 
lows: (i) The target nucleus is taken to be a collec- 
tion of independent nucleons distributed spatially 
according to a charge density profile p(r). (ii) The 
nucleons are regarded to be fixed within the nu- 
cleus neglecting Fermi motion and nucleon recoil 
effects. This implies that the pion energy E v does 
not change in the elastic scatterings. Furthermore 
the target is assumed to stay isotopically neutral 
(Z = N) during the multiple scattering processes. 

(iii) The pion interactions in the nucleus are as- 
sumed to be incoherent irN reactions taking place 
in the A resonance region. Since pion production 
and more complex channels are closed in this re- 
gion the rescattering processes can be described 
by two cross sections, the pion absorption cross 
section per nucleon a & hs{W) and the usual elas- 
tic wN cross sections. The relevant irN reactions 
read 

7T+ + N — > 7T+ + N, 7r + + n —>■ 7T° + p, 

7T° + N —> 7T° + N, 7T° + p -> 7T+ + 71, 
TT°+n^lT~+p, 7T~ + N -> 7T~ + N, 
7T~ + p — > 7T° + n, 

tt*' + N -> X, (tt £ X)(-> absorption) . 

(iv) In the A resonance region the irN cross sec- 
tion is dominated by the isospin 7 = 3/2 ampli- 
tude such that all elastic cross sections are related 
to ov+p(I / F) by Clebsch-Gordan coefficients. Let 

q% = (n i (7r + ),7i l (7r°),n i (7r~)) T (4) 

denote the initial multiplicity distribution of pi- 
ons in the medium (with fixed pion energy). In 
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a single elastic scattering this distribution will be 
modified to a final distribution qj\ 



<tf = Q Qi 



(5) 



with a known 3x3 matrix Q 1 which follows 
easily from a Clebsch-Gordan analysis of isospin. 

Accordingly, taking into account all possible 
multiple scatterings the final pion charge distri- 
bution is given by 



qf = Mqi 



M = J2 P nQ T 
n=0 



(6) 



where P„ is the probability that the pion exits 
the nucleus after exactly n irN scatterings. It 
should be noted that the matrix Q does not in- 
clude any absorption which is taken into account 
by X^^lo^™ < 1- The eigenvalues and eigenvec- 
tors of the matrix Q are given by 

Ax = l J g 1 = (l,l J l) T > (7) 
A 2 = f,g 2 = (l,0,-lf, (8) 
As - h, 93 = (1,-2, If (9) 



and one finds immediately the three eigenvalues 
of the charge exchange matrix M in dependence 
of the three eigenvalues A^ of Q 



/(A fc ) = ^P n A£,fe = 1,2,3 



(10) 



n=0 



The connection between the eigenbasis and the 
canonical basis is given by 

A p (l-c-d) = i/(l) + i/(|) + i/(i), 
A p d=i/(1)-I/(I), 
^V=i/(l)-i/(f) + i/(!) (11) 

or inversely 

/(I), 



Ap — 



|/(|)//(l) + s/(f )//(!), 



d = | [1 -/(!)//(!)] • 2 (12) 

As already mentioned in the introduction, the 
Pauli suppression factor of step 1 will be con- 
veniently absorbed into the normalization of the 
charge exchange matrix M 



A,, 



A p — . 



>(W,Q 2 )f(l) 



(13) 



The function /(A) contains the dynamical details 
of pion multiple scattering in the nucleus. An 
outline of how this function is calculated by solv- 
ing a transport problem for pions in the nucleus 
will be given in the next section. 

2.1. Calculation of /(A) 

To a very good approximation the three- 
dimensional problem can be reduced to a one- 
dimensional transport problem by projecting the 
forward-hemisphere of the scattered pion onto the 
forward direction and the backward-hemisphere 
onto the backward direction [see Eq. I|37jl]. In this 
approximation the pion scatters back and forth 
along its initial direction of motion until it is ab- 
sorbed or leaves the nucleus. As is explained in 
more detail in Refs. |1I2| the nucleon density pro- 
file along the scattering line can be scaled out of 
the problem by an appropriate change in length 
variable such that it is equivalent to consider a 
uniform one-dimensional nuclear medium extend- 
ing from x — to x — L with effective length 
(optical thickness) 



L = L{b) = 



1 



P(0) 



+ OC 



dz p(x/z 2 + b 2 ) , (14) 



where b is the impact parameter. Furthermore, 
to simplify the discussion we only consider for- 
ward scattering of the pions (and also neglect 
Pauli suppression in this step 2). The solutions 
taking into account forward and backward scat- 
tering, f± (A) , can be found in Appendix A of 1 . 
Apart from being more realistic, they are impor- 
tant for estimates of the background to proton 
decay by neutrino production of pions in the A 
resonance region 

The transport process can be described in 
terms of the basic probabilities for a pion density 
to propagate from y to x and to a) interact (scat- 
tering or absorption) b) scatter c) be absorbed in 
[x, x + dx] 



a) (x\ Ptot | y) 

b) (x\ P ccx \y) 

c) {x\ Pabs \y) 



Ke-w-'iQfy - x) (15) 
fi {x\ P tot \y) (16) 



a {x\ P to t | y) 



(17) 



with k = ;o(0)<7tot ('inverse interaction length'), 
/i = cr ccx /<7tot being the probability that the pion 
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is scattered and a = cr a bs/ct t the probability 
that the pion is absorbed (in a single scattering 
process). The choice of p(0) as density of the uni- 
form one-dimensional medium (i.e. k = p(0)at o t) 
is related to the normalization of the effective 
length in Eq. l(Tljl. Note also that fj, + a = 1 
since ctot = cr ccx + er a b s - The cross sections for 
charge exchange (cex) and absorption (abs) will 
be specified below. 

The density of pions in the medium after n 



scatterings, \pl 

sit y IPin } ) b y 



(n) 



is related to the initial den- 



» 



P 



cex rm 



(18) 



where the initial density is normalized to 
( x \P <: m'^ = 1/L such that the number of pions 
in the medium is one: 

JVff^jf d* (x\$) = l. (19) 
The density of pions leaving the medium after n 



scatterings, is given by 

(i-^tot) \ P [: ] ) 

(1 - Ptot)P< 



, (n) 
Pout 



J 0) ) 

cex I fin / 



(20) 
(21) 



The probability that the pion leaves the medium 
after exactly n rescatterings is then 

r-L 



Pn 



N. 



(n) 



da; {x\p^ t 



(22) 



The dynamical function /(A) we wish to calculate 
is then formally given by 



/(A) 



]Tp„a* 

n=0 



with 



n=0 



dx (x\ip out ) (23) 



(24) 



Using 1 - P tot = 1 - i + i ( 1 - APcex) with a:=Xp 
and £~ =0 A n P c ^ x = (1 - AP^)- 1 =: 1 + F we 
can write IV'out) as 

|Vout>= [1 + (1-|)F] \$>) . (25) 



Inserting Eq. mt o Eq. we find 

/(A) = l + (l-i)^ dx dy <a:| F |y> ± (26) 
where we have inserted the unit operator 

1= [ dy \y)(y\ (27) 



(i 



and used = \/L. 

Using the definition of the operator 1 + F — 
(1 - APcex)- 1 and (1 - APcex)- 1 ^ - AP CCX ) = 1 
we can furthermore write 



F — APcox + F AP C( 



(28) 



which leads us, using Eqs. 1)15(1 and 1(16(1 . to the 

transport integral equation 

f(y) = a(l-e- Ky )+(XK f dz f(z)e- K{y - z) (29) 
Jo 

with 



dx (x\ F \y) 



(30) 



The integral equation can be transformed into a 



differential equation by evaluating 4-[e KV f(y)\ for 



the left and the right side of Eq. ((29(1 resulting in 
d 



f= K (a-l)f(y)+Ka,f(0)=0. (31) 



This differential equation is solved by 
f{y) = ^—[l-h{y)/hm , 

1 — a 



(32) 



where h(y) — e K ^ CT_1 ^ is a solution of the ho- 
mogeneous equation h! = k(<t — l)h. The final 
solution for the dynamical function /(A) is now 
easily found to be 

/(A) = l + dy/(y) 

1 _ g-Ki(l-er) 

= kL{1-o) ' g = A/i - (33) 

Finally, we average the solution in Eq. 1(33(1 over 
impact parameters: 



/(A) = 



/ °° b db L(b)f(X, L(b)) 
f™bdbL(b) 



(34) 
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0.25 



0.225 - 




r(fm) 

Figure 1. Charge density profiles for oxygen, ar- 
gon, and iron normalized to J d 3 rp(r) = A. 



2.2. ANP model: Input 

In this subsection we summarize the input 
quantities entering the ANP model. 

(i) First of all the nucleon density p(r) enters 
the calculation of the effective length L(b) via 
L(b) = J dz p[r(z, &)]. The charge density 
profiles have been determined in electron scatter- 
ing experiments. For lighter targets like sO 16 the 
charge density is parameterized by a 'harmonic 
oscillator form' 

2 4 

p(r)=p(0)exp(-r 2 /i? 2 )(l+C^+C 1 ^ i ) (35) 

whereas for heavier targets like isAr i0 and 26-Fe 56 
a 'two parameters Fermi model' is utilized: 

p(r) = p(0) [1 + exp ((r - _1 . (36) 

The parameters p(0), R, C, and C\ can be found 
in Table 1 of 0. The corresponding charge den- 
sity profiles p(r) are shown in Fig. ^ 

(ii) The second input is the usual cross section 
for elastic ttN scattering in the (3, 3) resonance 
region. This region is dominated by the I = § 
amplitude. The elastic cross sections per nucleon 
for a ir k (k = ±, 0) in the initial state are given 
by Z/A (Tjckp + N/A cr^kn. For an isoscalar tar- 
get [Z = N) and neglecting the I = | channel 
they are independent of the pion charge and will 
be denoted by er ccx . Moreover, the elastic cross 
section cr cox is proportional to the cross section 



^ cx (x w+p (W) (l + 3cos 2 0) , (37) 

a™* = fow+ P (W) (38) 

with 

a« +p (W) = <r (3 , 3 ) (W) + 20 mb (39) 

where 17(3,3) (W 7 ) is a resonant cross section which 
can be found in £Q and the second term is a 
constant non-resonant background. In Eqs. I|37|l 
and l|38|l we have omitted for simplicity a possi- 
ble Pauli suppression factor taking into account 
the Pauli exclusion principle in step 2. Note also, 
that it is debatable whether such a factor should 
be included here since we are working in a picture 
with fixed nucleons. 

(iii) The third ingredient is the cross section per 
nucleon for pion absorption, a a bs{W). For this 
quantity Sternheim and Silbar have published two 
parametrizations 9 10 which they have extracted 
from data on single pion production in pA scat- 
tering. These parametrizations can also be found 
in pQ . In Fig. [5] these two absorption cross sec- 
tions, model (A) and model (B), are shown to- 
gether with the elastic and the corresponding to- 
tal cross sections. As can be seen the models (A) 
and (B) for <7 a b s are quite different in shape and 
normalization. 

3. Fixing the normalization of a a b s (W) 

The two models for the absorption cross sec- 
tion, model (A) [5j and model (B) ^Uj result in 
quite different total amounts of absorbed pions: 
for example, model (A) predicts that about 19% 
of the pions will be absorbed in oxygen compared 
to 43% in model (B). 

These numbers have been obtained in an 'av- 
eraging approximation' pQ in which the In- 
dependence of /(A, W) is averaged over the (3, 3) 
region. In this case the dynamical functions (the 
charge exchange matrix) are mainly sensitive to 
the region around W ~ toa and thus mainly sen- 
sitive to the normalization of cr^^W ~ wa). 

In the following we fix this normalization by 
using experimental data for 1-pion production in 
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Figure 2. Total, elastic, and absorption cross sec- 
tions in the (3,3) region. 



CC z^-deuteron [<-> free case] and ^ M -neon scat- 
tering. These data have been weighted to the 
same atmospheric E v spectrum in a paper by 
Merenyi et al. JI] allowing for direct comparison. 
In Table 1 of 11 relative charged current popula- 
tions (w.r.t. the total cross section) are provided 
for deuteron (D) and neon (Ne). The fractional 
contributions of the 1-pion production channels 
(7r + , 7T°, 7r~) are given by: 



D : n 
Ne : ff 



(0.165, 0.09, 0) T (40) 

(0.11, 0.05, 0.01) T 

±(0.014, 0.02, 0.01) T . (41) 



The fractional populations fi and ff are related 
in the ANP model by: 

f,=mxK, K= 1 -^§. (42) 

Unfortunately the total cross sections, weighted 
to an atmospheric neutrino flux, have not been 
published in 11 . However, the correction factor 
K should be close to one since the total cross 
section for Ne is not affected by charge exchange 
and pion absorption. In the following exercise we 
assume K = 1 and later assign a normalization 
uncertainty of 3%. 

Furthermore, it is not viable to solve ff = 
M[A P , d, c] fi for the three parameters A p ,d, c 



since the solution strongly varies within the er- 
rors of ff. Instead, we have made the reason- 
able assumption < c < d which means that 
the probability for it" — > n + is smaller than 
the probability for 7r° — > ir + . Under this as- 
sumption one can fit the two parameters A p 
and d for a fixed value of c as long as c < d: 
c A p d 

c = ^FiE 0695 0.147 

c = 0.01 -> Fit: 0.696 0.128 

c = 0.02 -> Fit: 0.696 0.109 

c = 0.03 -> Fit: 0.696 0.091 

c = 0.04 -> Fit: 0.697 0.072 

c = 0.05 -> Fit: 0.698 0.053 
[x 2 /d.o.f ~ 0.4] 

As can be seen the parameter A p is well con- 
strained: A p — 0.696±0.002. More conservatively 
we use in the following A p = 0.70 ± 0.02 taking 
into account the normalization uncertainty due 
to the correction factor K . On the other hand 
the parameters d and c are correlated and adopt 
values in the range c G [0,0.05], d G [0.15,0.05]. 

Inspecting Eq. (fTTTf) we find that /(A = 1) = 
12^=0 Pn = ^ — A where A is the fraction of ab- 
sorbed pions. Therefore, the averaged fraction A 
of absorbed pions is given hy A — 1 — /(A = 1) 
where /(A) is the averaged dynamical function 
which is related to the above extracted parame- 
ter Ap via Ap = g(W, Q 2 )f{l). In order to es- 
timate the effect of the Pauli suppression factor 
we take W ~ ttia and Q 2 ~ 0.1 GeV 2 leading to 
g(W,Q 2 ) = 0.93 ± 0.05 [see Table 3 in 0]. Us- 
ing these values we finally find /(l) = 0.75 ± 0.05 
implying that about 25% of the pions have been 
absorbed in the neon target: A = 0.25 ± 0.05. 

A fraction of 25% pion absorption can be ob- 
tained by renormalizing absorption model (B) 
[model (A)] by a factor ~ 0.3 [~ 1.4]. With this 
renormalization we can compare the ANP model 
with the above fitted results. The ANP model 
gives for both renormalized absorption models 
for a neon target /(l) = 0.75 (by construction), 
d = 0.15, and c = 0.05. This result compares 
favorably with the fitted values. 
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4. Dynamics of the ANP model 

4.1. Linearisation 

The simple solution in Eq. (|33[1 is well-suited 
for further analytical investigations of pion multi- 
ple scattering. For example, the forward solution 
jnHJl at A = 1. 



f(X = l,L,W) 



1 



PO La abs 



(43) 



can be considered in the limit poLo^s <C 1 relat- 
ing the fraction of absorbed pions A to the ab- 
sorption cross section per nucleon a a bs(W): 

A{L, W) = 1-/(1, L, W) ~ \p L(j &hs {W) . (44) 

Averaging L = L(b) over the impact param- 
eters b for oxygen as target material gives L ~ 
1.9R with radius R ~ 1.833 fm. The nuclear 
density for oxygen, po = 0.141 fm -3 , has been 
taken from Table 1 in Replacing pqL by 

poL ~ 0.05 mb -1 in Eq. (|44|) and taking into ac- 
count the averaging over the (3, 3) resonance by 
the replacement W — > tua we arrive at the fol- 
lowing rule of thumb for the fraction of absorbed 
pions in oxygen 



A ~ 0.025 a ahs (W ~ m A )[mb] . 



(45) 



The absorption cross sections can be taken 
from Fig. [2] and we find the following fractions 
of absorbed pions: 

model <7 a bs(W = m A ) A 

(A) ~ 6.0 mb 15%r~ 

(B) ~ 28.4 mb 71% 
0.3 x (B) ~ 8.5 mb 21% 

Here the renormalization factor 0.3 for model (B) 
has been taken from the preceding section. The 
same result is found by renormalizing model (A) 
with a factor 1.4. Note also that the 71% for (the 
original) model (B) is not realistic because the 
condition poLa^bs <C 1 is not satisfied. 

4.2. 'How many' multiple scatterings? 

Figure[3]shows the dynamical function f(X, W) 
for oxygen in dependence of A for several values 
of W. Writing /(A) = XT=o p n X n with the prob- 
abilities Pk (k — 0, 1, 2, . . .) that the pion is ob- 
served after k ttN scatterings it is obvious that 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
X 

Figure 3. Dynamical function /(A, W) in depen- 
dence of A for several values of W . 



f(X = 0) = P and /(A = 1) = 1 - A where A 
is the probability for pion absorption. [Note that 
Sfc°=o Pk + A = 1.] It is an interesting question, 
how large the probabilities P n , n > for multiple 
scattering are. Qualitatively, a stronger curvature 
of the function /(A) indicates a higher probabil- 
ity for multiple scattering. As can be seen, in 
the vicinity of W = toa the probability that the 
pion rescatters several time is largest whereas at 
W = 1.1 GeV the curve is almost linear such that 
only Po and Pi contribute appreciably. 

Of course the probabilities P„ can be calculated 
exactly within the ANP model according to Eq. 
H22fl or by differentiating the solution for /(A) in 

Eq. Pn = ^S/|A=0-^ 

On the other hand, inspired by the fact that 
the pion energy remains constant, it is interesting 
to make the assumption (which is only asymp- 
totically correct) that Pfc+i ~ rPj, for k > 1 
or equivalently Pk — Pit' 0-1 for k > 1. Under 
this assumption /(A) is a geometrical series and 
is given by the simple solution 



/(A) = ^P„A n cPo + T 



n=0 



PA 
-rX 



(46) 



The constant r can be fixed from /(A = 1) = 1— A 
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Table 1 



Probabilities in % for multiple scattering. 



WfGeVl 


1.1 


1.2 


1.25 


1.5 


A 


0.0 


11.8 


25.7 


48.4 


Po 


90.8 


33.6 


28.8 


43.4 


Pi 


8.0 


14.1 


12.2 


6.5 


P2 


1.0 


10.5 


8.9 


1.3 


p 3 


0.13 


7.8 


6.5 


0.26 


Pi 




5.8 


4.8 




p 5 




4.3 


3.5 




p 6 




3.2 


2.6 




Pj 




2.4 


1.9 




Ps 




1.8 


1.4 




P9 




1.3 


1.0 




P10 




1.0 


0.7 





resulting in r = 1 — 1 _ A 1 _ Pg ■ 

Thus, once A, Pq and Pi are known, all higher 
probabilities and the function /(A) can be easily 
estimated. The result of such a procedure is listed 
in the following table where Pq, Pi and A have 
been determined from Fig. [3J 
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5. Conclusions 

In this article, we have given a basic intro- 
duction to the pion multiple scattering model 
of Adler, Nussinov, and Paschos focusing on 
the input parameters of the model, particularly 
on the cross section for pion absorption cr s ^ s (W) 
which is the least well determined ingredient. Us- 
ing data for i/ M -deuteron and i/^-neon scattering 
|1 1 j we could fix the normalization of <r a bs(W0 
corresponding to a fraction of (25 ±5)% of ab- 
sorbed pions in neon. The parameters A p , d, and 
c have been determined to be A p = 0.70 ± 0.02, 
c £ [0,0.05], and d £ [0.15,0.05] which compares 
quite favorably with the predictions of the ANP 
model for neon. 

In order to test and improve the ANP model it 
will be necessary to make detailed measurements 
of single pion electroproduction in the region of 
the (3, 3) resonance using different heavy targets 
and to compare it with the corresponding cross 
sections on free nucleons 0. 



